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We classify all possible implementations of an abelian symmetry in the two-Higgs-doublet model 
with fermions. We identify those symmetries which are consistent non-vanishing quark masses and 
a CKM matrix which is not block-diagonal. Our analysis takes us from a plethora of possibilities 
down to 246 relevant cases, requiring only 34 distinct matrix forms. We show that applying Z n with 
n > 4 to the scalar sector leads to a continuous (7(1) symmetry in the whole Lagrangian. Finally, 
we address the possibilities of spontaneous CP violation and of natural suppression of the flavour 
changing neutral currents. We explain why our work is relevant even for non-abelian symmetries. 
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I. INTRODUCTION 



The least known aspect of the electroweak interactions is its scalar sector. In the Standard Model (SM) there is 
r— i only one Higgs but, although this is an economical choice, there is no fundamental reason for nature to adopt it. 
'■^ ' Ultimately, the number of Higgs fields, like the number of fermion families before it, must be assessed experimentally. 
""P" 1 . Partly for this reason, there has been a great interest in multi-Higgs models. This is also due to the fact that many 
O i' interesting new effects arise, such as the presence of charged scalars, the possibility for CP violation in the scalar 
sector, and the possibility for spontaneous CP violation, to name a few. 

One problem with multi-Higgs models is that they involve many more parameters than needed in the SM. This 
\ problem can be tamed by invoking discrete symmetries. A complete classification of the impact of discrete and 
t-H ■ continuous symmetries in the scalar sector of the two Higgs doublet model (THDM) has been discussed in the 
literature [l|, , and some incursions exist into theories with more than two Higgs doublets [!, H[ . There are also 
several articles discussing specific implementations of discrete symmetries in both the scalar and fermion sectors, but 
' no complete classification exists. This is the problem we tackle here. 

This article is organized as follows. In section [TT] we introduce our notation and show the impact that a choice of 
abelian symmetries in the scalar and fermion sectors has on the Yukawa matrices. A priori there are 3 18 possibilities. 
In section UHl we show how simple experimental considerations, such as the absence of massless quarks and the non- 
block-diagonal nature of the CKM matrix can be used to curtail this number down to 246. Up to permutations, these 
involve only 34 forms of Yukawa matrices, which we show explicitly. Since any finite discrete group has an abelian 
sub-group, our classification is important even for those considering non-abelian family symmetries. We present two 
important results in section IIV1 Our classification is then used to address two questions: whether one can have 
spontaneous CP violation, in section [V] and whether one can relate the flavour changing neutral current interactions 
5_i ■ with the CKM matrix, in section IVI1 We draw our conclusions in section IVIII 

a _ 

II. NOTATION 
A. The Lagrangian 

Let us consider a SU{2) ® U(l) gauge theory with two hypercharge-one Higgs-doublets, denoted by $ a , where 
a = 1,2. The scalar potential may be written as 

- £ H = = Y ab (^ b ) + iZ o6 , cd ($t$ 6 )($t$ d)i (1 ) 

where Hermiticity implies 

Y a b = Y ba , 

Z a b,cd = Zd^ab = Z badc . (2) 

Minimization of this potential leads to the vacuum expectation values (vevs) ($ a ) = v a . 
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The theory contains also 3 families of left-handed quark doublets (ql), right-handed down-type quarks (tir), and 
right-handed up-type quarks (pr). For the most part, we will ignore the leptonic sector, since the analysis would be 
similar. The Yukawa Lagrangian may be written as 



Cy = -q L 



(ri$i + r 2 $ 2 W + (Ai$i + A 2 ^Pr +H.c, (3) 



where = iT2^* k , and qL, nR, and pr are 3- vectors in flavour space. The 3x3 matrices F^, A&, contain the complex 
Yukawa couplings to the right-handed down-type quarks and up-type quarks, respectively. 

B. Basis transformations 

The Lagrangian can be rewritten in terms of new fields obtained from the original ones by simple basis transfor- 
mations 

$ Q -> K = u ab $ 6 , 
ql -> q' L = Ul qL, 

n R ~> n R = U n R. n R , 

PR -> PR = UpR PR, (4) 

where U £ U (2) is a 2 x 2 unitary matrix, while {Ul, U„r, U p r} € U(3) are 3x3 unitary matrices. Under these unitary 
basis transformations, the gauge-kinetic terms are unchanged, but the coefficients Y ao and Z ao ^ ca \ are transformed as 

Yab Y ab = U aa Y a p U b p, (5) 

Zab.cd — > Z' abcd = U aa U cl Z a p n s UfrpUjg, (6) 

while the Yukawa matrices change as 

r a ^ r; = u L r a u\ R (c/t) Qa 

A Q A' a — Ul Aq, Up R (U T ) aa . (7) 

Notice that we have kept the notation of showing explicitly the indices in scalar-space, while using matrix formulation 
for the quark flavour spaces. The basis transformations may be utilized in order to absorb some of the degrees of 
freedom of Y, Z, T, and/or A, which implies that not all parameters in the Lagrangian have physical significance. 

C. Symmetries in the THDM 

We will now assume that the Lagrangian is invariant under the symmetry 

* a -> $f = S ab $b, 
qL -> q'[ = S L qL, 
nR -> n s R = S nR n R , 

PR -> Pr = S pR PR-> ( 8 ) 
where S € U{2), while {Sl, S n R, S p r} £ U(3). As a result of this symmetry, 

Y ao — S aa Y a p S b p, (9) 

Zab,cd — Saa Scy Z a fj_^s Sfrp S^g, (10) 

F a = S L T a Sl R (5t) Qa , (11) 

A a = S L A a Sl R (S T ) aa . (12) 

Under the basis transformation of Eq. (0} , the specific form of the symmetry in Eq. ([5} is altered as 

S' = U S U\ (13) 

S' L = U L S l u[, (14) 

S'nR = U n R S n R U nR , (15) 
S'pR = U p R S p R U pR , (16) 
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Suppose that one has chosen to apply the symmetry {S, Sl, S nR , S pR } in some basis. By a judicious choice of 
{£/, Ul, U nR , U pR } one may bring the symmetry into the form 

S = diagje^e^}, (17) 

Sl = diag{e iQ1 ,e lQ2 ,e ia3 } , (18) 

S nR = diag{e^,e^,e^}, (19) 

S pR = diag{e^,e^,e^}. (20) 

What about global phases? Clearly, an overall phase change has no effect on the symmetry. For example, taking 
U = e t6 \i, leaves S' = S. However, it is easy to see from Eqs. (|9"1)- (fT2"|) that the symmetry 

S = e iS S, S L = e i& S L , S nR = e^S nR , S pR = e^S pR , (21) 
imposes the same restrictions on the Lagrangian as the symmetry {S, Sl, S nR , S pR }, as long as 

e i0-&-§) = 1 and e i(7-5+e") = i (22) 

This can be used to bring Eqs. (fTT)) - (f2"U f into the form 

S = diag{l,e i9 }, (23) 

S L = diag{e iQ1 ,e iQa ,e iQ3 } , with a x = 0, (24) 

S nR = diag{ e ^,e*,e 4fe }, (25) 

S pR = diag{e^,e^,e^}. (26) 

For 6 — it, S — diag(l, —1) leads to the usual Z2 Higgs potential. Any other value of < 9 < 2ir, leads to the full 
U(l) symmetric Higgs potential. For example, with 9 = 2ir/3, S 3 — I2, and a Z$ symmetry is imposed on the scalar 
fields. Nevertheless, because the scalar potential only has quadratic and quartic terms, the resulting Higgs potential 
has the full U(l) Peccei-Quinn symmetry Q. If this symmetry is broken spontaneously by the vacuum, we will have 
massless particles. As a result, great care must be taken when imposing what may look like discrete symmetries in 
multi-Higgs models. Substituting Eqs. (|T7 l) -([2"0 ]l in Eqs. {TTJ and ([11]), we find 

(r B )« = e^-ft-^(r a ) ii; (27) 
(A ) y = e^'-^HA^, (28) 

where no sum over i and j is intended on the right-hand sides. For the simplified form in Eq. (|23[) we set 6\ = 
and 92 = 9. Furthermore, we will always take 9 7^ (mod 2ir), since we are only interested in symmetries which do 
transform the scalar fields. It will prove useful to keep a\ explicitly, bearing in mind that it can be set equal to zero 
without loss of generality. These equations constitute our starting point for what follows. 



D. Preliminary constraints on the Yukawa matrices 

We will concentrate first on the down-type Yukawa matrices T a . Given a symmetry written in the form of Eqs. (|23[) - 
([2"6]) we conclude from Eq. d^ZJ) that 

• (ri)^ can take any value if fty = 0; 
. (r^. =0if^0; 

• (-^2)^ can take any value if 9ij = 9; 
. (T 2 ) tJ = if % ^ 9; 

where we have defined 



9ij = a» - /3j. 

We conclude that, for a matrix S characterized by a given 9 ^ 0, there are only three possibilities: 



(29) 
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1. 







(Ti 



any 



2. 9ij — 



(r 2 



o 



(rx) y = and (T a ) tf = any; 

3.0^0,0 => (rxjy = = (r 3 ) y . 

All conditions on 0y are mod(27r). Noticing that only five 
6*2i, and 03 1. Then, 



are independent, we will take these to be On, 6*12, 6*13, 



?22 
732 



^21 

hi 



'12 

?12 



'11 j 



723 

733 



721 
931 



?13 
?13 



711, 
711- 



(30) 



For each 7^ 0, we must only consider five 9{j. The possibilities #y = and 0y = are simple to enumerate. 
Unfortunately, the impact of 0y 7^ 0, 9 depends on the exact value of 6ij. Thus, there are far more than the 3 5 
possibilities one might naively expect. For example, choosing {0n, 612, #21} = {79, 29, 29} and 9 — V2tt, we conclude 
that the (1,1), (1,2), and (2,1) entries of Ti and T2 matrices vanish, as do the (2,2) entries. In contrast, choosing 
{0n, 612, 621} = {40, 20, 20} and = V2n we conclude that the (1,1), (1, 2), and (2, 1) entries of Ti and T 2 matrices 
vanish, but the (2, 2) entry of Ti need not vanish 1 . 

Some possibilities are trivially inconsistent with experiment. For example, choosing {0n, 0i2, 0i3, 021, 03i} = 
{0, 0, 0, 0, 0}, then the matrix 



becomes 



e 



e = {Bn} 



9 9 
9 29 29 
9 29 29 



(31) 



(32) 



For 0^0 (identity operation) and 9 7^ 7r (usual Z 2 symmetry), we are lead to Yukawa matrices of the form 



Ti = 



an 





" 







bi2 


b 13 











, r 2 = 


&21 



















631 









(33) 



Upon spontaneous electroweak symmetry breaking, the down-type quark mass matrix will arise from the bi- 
diagonalization of 



V1T1 + v 2 T 2 = 



viau v 2 b 12 v 2 b 13 
W2&21 
v 2 b 3 i 



(34) 








I>i2 


013 




an 





" 


ri = 


021 








, r 2 = 













_ &31 




















whose determinant is zero. As a result, this model would lead to one massless quark, which is ruled out by experiment. 
Notice that choosing {0n, 0i2, 013, 021, 03i} = {0, 0, 0, 0, 0} would lead to Yukawa matrices of the form 



(35) 



This is the same as Eq. (1331) . with the substitution $1 o $2- Said otherwise, these possibilities represent the same 
model. The interchange $1 <-> $2 cuts down the number of distinct models by almost a factor of two. 
An old model by Lavoura [5J| had 

S = diag{l,-l}, 
S L = diag {1,1,1}, 
S nR = diag {1,1,-1}, 

S pR = diag {1,1,1}. (36) 



1 Notice that the freedom to choose 01 = does not reduce the number of possibilities. 
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Thus 



e 



o o o 
o o e 
o o e 



(37) 



leading to Yukawa matrices of the form 





an 


a\2 


o" 




" 





bi3 


ri = 




a-22 





, r 2 = 








&23 




_ a 3 i 


032 













^33 _ 



(38) 



A model where 



e = 



6*0 
6*0 
6*0 



(39) 



will be indistinguishable from Lavoura's model, as will a model where the 0s move to the first column. Such permu- 
tations will further cut down the number of distinct models. 
For the up-type sector we define 



on 



(40) 



As before, for a matrix S characterized by a given 7^ 0, there are only three possibilities: 
=> (Ai) tf = any and (A 2 ) ij = ; 
=> (Ai) tf =0 and (A 2 ) ij = any; 
=► (A 1 ).. = = (A 2 )... 



1. 9 l3 = 



2. 9ij — —9 



3. Sij + 0, 



All conditions on 9^ are mod(27r). Clearly we can choose independently 9u, 0i 2 , and #13, and then 



721 

hi 



721 
931 



011 + 011 
#11 + 011 



022 
032 



021 
031 



011 + 012, 
011 + 012, 



?23 
^33 



721 
731 



011 + 013, 
011 + 013- 



(41) 



There are 9 entries in the down-type Yukawa matrices. For each there are only three possibilities (the entry exists in 
Ti but not in T 2 ; the entry exists in T 2 but not in Ti; the entry does not exist in either). The same occurs for the 
up-type Yukawa matrices. As a result, we would have potentially 3 18 possibilities. But, as we have illustrated above, 
interchange and permutations help cut this number down. More importantly, many of the models entail massless 
quarks, a diagonal CKM matrix, or other inconsistencies with experiment. These are ruled out. This is what we turn 
to next. 



III. MODEL CLASSIFICATION 



A. The left-space 

The left-handed space (where the left-handed quark doublets live) is rather constrained because it affects the down- 
type quark mass matrix, the up-type quark mass matrix, and also the CKM matrix. The quark mass matrices are 
obtained by bi-diagonalizing the matrices 

r ee «iri + «2r 2) (42) 

A ee t£Ai+t£A 2 ,. (43) 
whose two indices live in different spaces. But both indices of the hermitian matrices 

H d ee rrt = |«i| 2 rir{ + M 2 r 2 r 2 + v^r^l + ^ 2 r 2 rl (44) 

H u ee AA f = |w 1 | 2 A 1 A t 1 + |i; 2 | 2 A 2 A 2 +i;*w 2 A 1 A^ + w 1 i> 2 *A 2 Al (45) 
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live on the left-space. These matrices can be diagonalized through unitary matrices VdL and V u l as 

V dL H d VJ L =D 2 d = diag{m 2 d , ml, m 2 b }, 
V uL H u vl L = D 2 U = diag{m 2 u ,m 2 c ,m 2 t }, 



(46) 
(47) 



where V = V uL Vj L is the CKM matrix. 

We may now see the impact of the symmetry on the left-space and how it affects the quark masses and mixings. 
We start from Eq. (jlip in the form 



Tl — Sl Tl SnR, 

T2 = Sl ^2 Sl 



J nR 



which, using the simplified form of Sl in Eq. (1181) . we can combine into 

r x rl = s L r x rj s[ = 



An 
A 2 i e~ 



A 12 e Ml2 

122 



A 3 i e^ 1 A 32 e- 



A 13 e 

^33 



r 2 r| 



Sl ^2^2 &L 



T.Tl = S L riTl S{ e l6 = 



v 2 v\ = s L t 2 t\ sl e 



B u B 12 e Ml2 Bis e~ M31 
B 21 e" M12 B 22 B 23 e JQ23 
Bai B 32 B 33 



C n e xd C u e l{a ^-+^ Ci 3 e^" 31 ^ 
C 21 e - i( - ai2 -^ C 22 e ie C 23 e l ^ +e ^ 
C 31 e l ( Q31+0 ) C 32 e - 1 ^^ C 33 e l6 



D 12 e 



i(ai2—t 
40 



D l3 e 



-i(a 3 i+9) 



D 21 e -^+ s ) D 22 e 
D 31 e l< - a31 -^ D 32 e-^^+e) 



D 23 e 
D 33 e 



i(a 2 3-S) 



(48) 



(49) 



(50) 



(51) 



(52) 



In the previous four equations, A = TiT\, B = T 2 V\, C = TiY\, and D — T 2 V\, respectively. We have defined 

ai2 = ai — a 2 , a 23 = a 2 ~a 3l a 31 = a 3 - at, (53) 

which satisfy 

ai 2 + a 23 + a 3 i = 0. (54) 



It is easy to see that the up-type Yukawa matrices satisfy identical equations, with 9 — > — 9. 
We define the set 

J = {x : x = (mod 2tt) V x = 9 (mod 2tt) V X = -9 (mod 2tt)}. 



(55) 



If ai2, a 23 , a 3 i £ J, then the matrices TiT\, T 2 t\, AiA{, and A 2 Aj are diagonal, while all 12 and 21 combinations 
vanish. As a result, Hd and H u are diagonal and the CKM matrix V is unity. This is ruled out by experiment. As 
a result, at least one must belong to J . Let us imagine that a\ 2 € J, while 023,031 ^ J. In that case, Hd and 
H u are block diagonal, so are the matrices VdL and V u l, and so is the CKM matrix V. This is also ruled out by 
experiment. We are left with the cases where 

1. one Qy is not in J ', while the two others are in J; 

2. all cky belong to J . 



Next we study these cases in detail. 
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B. Odd one out 



We look at the case where only one Qy is not in J . Let us take this to be «i2 £ 3 ', ct23, «3i <E 3 ■ It is easy to see that 
the only possibilities that satisfy this requirement are {«i2, Q!23, 031} = {2#, —9, —9} and {ai2, 023, «3i} = {—29, 9, 9}. 
The second possibility arises from the first through the interchange ol\ 002. These symmetries act on the left of the 
Yukawa matrices and, thus, we go from one to the other by simply interchanging the first two rows of the corresponding 
Yukawa matrices. Similarly, the relevant cases where «23 ^ 3 ' 1 o^ii^is € 3 ", and a^i (jz. 3 , 0^2,0^23 6 J are related 
to the case shown here by mere permutations among the rows of the respective Yukawa matrices. As a result, we 
show only the case {«i2, "23, "31} = {2#, —#,—#}• Using Eqs. (|2"5j) . we obtain 9 2 i = #11 — 29, 6 3 i = On — 9. From 
Eqs. flU) we get 



e 



(711 (712 Cl3 

7ii — 29 9\2 — 29 $13 — 29 
#11 — 9 9\2 — 9 $13 — 9 



(56) 



The entries of this matrix which equal (mod 2-7r) lead to corresponding entries in Li; those which equal 9 (mod 27r) 
lead to corresponding entries in F2; all others lead to vanishing entries in Fi, T2, and, thus, in T. Recall that T 
cannot have a row of zeros nor a column of zeros; otherwise there would be a massless quark. This is a very powerful 
constraint. Let us consider the columns first. Since there must be at least one entry on each column, we conclude 
that 9ij g {0, 9, 29, 39} (mod 2ir). This would seem to lead to 4 3 possibilities. However, if On = #12 = #13, then there 
would be a (forbidden) row of zeros. The reason for this is that we are considering the case where 29 = ai 2 ^ 3 ', 
implying that 9 ^ zitt and 9 ^ z 2 27r/3 with z\ and Z2 integers - keeping the interval [0,27r[, 9 £ {0, 27r/3, n, Att/3}. 
This means that —9, ±29, and 36* can never equal (mod 27r), nor can they equal 9 (mod 2tt). Consider for example 
the possibility that On = 9 i2 = #13 = 39. Then, 9 would have 39 on the first row, 9 on the second row, and 29 on 
the last row. Because 39 and 29 cannot equal nor 9 (mod 27r), this would imply that the first and last rows of Fi, 
F2, and r vanish, leading to massless quarks. Also, possibilities where two Oij are equal to or to 39 lead to a 2 x 2 
block of zeros in T (implying massless quarks) and are, thus, excluded. There remain only eight independent forms 
for the Fi matrices (9 £ {0, 2tt/3, tt, 4tt/3}): 



• #11 — #12 — 9 


#13 = 2# 


(mod 27r) 












Ti = 


X 

X X 


, r 2 = 


X X 

X 


, # ^ 27r/3,7r,47r/3; 2 


• #11 = #12 = 9 


#13 = 39 


(mod 27r) 












Ti = 


X X 


, r 2 = 


X X 

X 


, 9 £ 271-/3,71-, 4tt/3; 



(57) 



(58) 



• #11 = #12 = 29, #13 = (mod 2tt) 





X 






Fi = 


X X 


, r 2 = 










X X 



# ^ 27r/3,7r,47r/3; 



(59) 



?ii = #12 = 29, #13 = # (mod 2tt) 









X 


ri = 


X X 


, r 2 = 






X 




X X 



# ^ 27r/3,7r,47T/3; 



(60) 



2 Eqs. H57I I— (164 ft are invariant under the symmetries for all 8, but they are only the most general forms consistent with the symmetry for 
those symmetries where 9 27r/3, 7T, 47r/3. See section \l V B I for details. 
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• U = 0, 6»i2 = 20, 6»i3 = (mod 2tt) 





X 




X 


ri = 


X 


, r 2 = 






X 




X 



9 ^ 27r/3,7r, 4tt/3; 



(61) 



hi = 0, 6»i2 = 39, 0i3 = (mod 2tt) 





X 




X 


ri = 




, r 2 = 


X 




X 







^ 27r/3,7r,47r/3; 



(62) 



• 0ii = 0, 0i2 = 20, 0i 3 = 30 (mod 2tt) 





X 








r\ = 


X 


, r 2 = 


X 


, ^ 27r/3,7r,4vr/3; 








X 





(63) 



On = 0, 0u = 20, 0i3 = 30 (mod 2tt) 









X 


ri - 


X 


, r 2 = 


X 




X 




X 



^ 27r/3,7r,47r/3. 



(64) 



The x denote an allowed complex entry; vacant positions mean that the entry is zero. All other allowed cases with 
{ai2, «23, 031} = {20, —0, —0} are related to these by permutations among the columns. This corresponds to a mere 
renaming ol the down-type right-handed fields {rim, riRi, nm}, having no physical significance. As explained above, 
all permutations of the rows correspond to physically allowed cases other than {an, 023, 031} = {20,-0,-0}. As a 
result, all column and row permutations of the matrices in Eqs. (|57[) - (I64I) correspond to physically allowed models; 
permutations on columns have no physical effect; permutations on rows also have no physical effect but must be 
performed simultaneously on the down- type matrices T and on the up- type matrices A. 



C. All in J 



We now turn to the cases where 012,023,031 G J ■ This means that each o^ can only take the values 0, 0, or — 
(mod 2ir). There would seem to be 3 3 possibilities. But Eq. (I5H) allows us to exclude a few. For example, taking 
—&12 = «23 = 0:31 = (mod 27r) into Eq. (|54[) would mean that = (mod 2n), a case we are not considering 
since it corresponds to unconstrained scalar fields: $1 — > $1, <J> 2 — > $ 2 . There are some cases which are possible only 
for specific values of 0. Postponing those for subsections MI C 41 and MI C 51 we are left with the following cases: i) 
{ai2, a 2 3, «3i} = {0, 0, 0}; ii) {«i2, 023, 031} = {0, —9, 9} (interchanging rows on the Yukawa matrices for this case one 
reaches the cases {012,023,031} = {—9,6,0}) and {012,023,031} = {0,0,-0}); and iii) {012,023,031} = {0,0, -9} 
(interchanging rows on the Yukawa matrices for this case one reaches the cases {ai2, «23, 031} = {0,-9,0}) and 
{0:12, 0:23,031} = {-9,0,9}). 



1. {<*i2, c*23, «3i} = {0,0,0} and any 6 



In this case, a% — a-i — 03 and 0n = 02 1 = 031) leading to 



e 



11 


1712 


fl3 


11 


012 


013 


11 


012 


013 



(65) 



Because a column of zeros in both Ti and T 2 would lead to massless quarks, we must have 0ij G {0,0}. There are 
2 3 possibilities; each column must exist in either Ti or IV Ignoring cases which differ only by permutation of the 
columns, we are left with the following structures: 



All 0i j equal 



Ti = 



X 


X 


X 






X 


X 


X 


, r 2 = 




X 


X 


X 







any 0; 



• Two 0\j equal 



Ti = 



• One 9 1 j equals 



X 


X 






X 


X 


X 


, r 2 = 




X 


X 


X 






X 


X 






X 


X 


X 




, r 2 = 


X 


X 


X 






X 


X 



any 0; 



Ti = 

This is the same as Eq. (|67|) . with the interchange <£>i <-> $2. 
• No 9ij equals 



any ti. 







X 


X 


X 




, r 2 = 


X 


X 


X 






X 


X 


X 



This is the same as Eq. (|66|) . with the interchange $1 <-> $ 2 . 

2. {ai2 , «23 ,0131} = {0, — 8, 8} and any 

Here 3 9 2 i = On, 9 31 = n + 9, and 



any ti. 



e = 



Cll 1712 <713 

^11 012 013 

711 + 9l2 + 013 + 



(66) 



(67) 



(68) 



(69) 



(70) 



implying that 0\j € {0,9,-0}. Ignoring cases which differ only by permutation of the columns, we are left with the 
following structures: 



{011,012,013} = {0,0,0} 



{011,012,013} = {0,0,0} 



r : = 



X 


X 


X 










X 


X 


X 


, r 2 = 
















X 


X 


X 


X 


X 










X 


X 


X 




, r 2 = 






X 










X 


X 




X 


X 










X 


X 


X 




, r 2 = 






X 






X 




X 


X 





any ti; 



(71) 



= 7T. 



(72) 



(73) 



The cases with {9n, 6*i 2 , $13} equal to {0,0,0} and {0,0,0} are obtained from these through column permuta- 
tions. 



3 Recall that the cases {012,023,031} = {— 9,9,0}) and {012,023,031} = {9,0,-9} are obtained from this through permutations on 
the rows of the Yukawa matrices. 
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{011,012,013} = {0,0,9} 





X 






X 


X 


1 l = 


X 


1 




X 


X 








X 








X 






X 


X 


r x = 


X 


, r 2 = 




X 


X 




X X 




X 







e = n. 



(74) 



(75) 



The cases with {0n, 0i2, 013} equal to {8,9,0} and {0,0,0} are obtained from these through column permuta- 
tions. 

{011,012,013} = {0,0,-0} 





X 


X 






r x = 


X 


X 


, r 2 = 








X 




X X 



(76) 



Setting 9 — ir we re-obtain Eq. ([73]) . The cases with {0n, 012, 013} equal to {0, —0, 0} and {—0, 0, 0} are obtained 
from these through column permutations. 



{011,012,013} = {0,0,-0} 



Ti 



X 






X 


X 




, r 2 = 


X 




X 




X 



(77) 



Setting = 7T we re-obtain Eq. l[75]). The cases with {0 U , 9 12 , 9i 3 } equal to {0,-0,0}, {0,-0,0}, {0,0,-0}, 
{—0,0,0}, and {—0,0,0} are obtained from these through column permutations. 



{0n, 012, 013} = {0, 0, —0} 











X X 






r x = 






, r 2 = 


X X 




, O^n; 






X 


















X X 


X 




Ti = 






, r 2 = 


X X 


X 


, = 7T. 




X X 


X 










equal to {0, — 


9,9}, and {-6 


,0,0} 


are obtained from 



(78) 



(79) 



mutations. 

For those wishing to check that all possibilities have been considered, we refer to the footnote 4 . 



4 We have also checked that 

• The cases where {0n , 012 , $13} equal {0,-8,-6}, {—9,0,-8}, and {—8,-8,0} lead to vanishing quark masses, if 8 7^ it, and to 
Eq. (75J, if 8 = tt; 

• The cases where {#11, $12, 613} equal {9, 9, 9} lead to vanishing quark masses, if 9 ^ tt, and to Eq. (1796 . if 9 = tt; 

• The cases where {0n, 012, $13} equal {9,-9,-9}, {—9,9,-8}, and {—8,-8,8} lead to vanishing quark masses, if 9 ^ n, and to 
Eq. (22), if 9 = vr; 

• The cases where {0n, $12, $13} equal {—9, —8, —9} lead to vanishing quark masses, if 9 7^ it, and to Eq. I I79I I. if 9 = tt. 
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{ai2 


, «23, Q3l} 


= {0,0 


—8} and c 






012 


013 


e = 


011 


012 


013 




011-0 


012- t 


) e 13 -e 



Here 5 9 21 = 0n, 631 = &n - 0, and 



implying that 0y G {0, 0, 20}. 

Ignoring cases which differ only by permutation of the columns, we are left with the following structures: 

• {011,012,013} = {0,0,0} 



(80) 





X 


X 




X 


ri = 


X 


X 


, r 2 = 


X 






X 







(81) 



Performing $1 O $2 and exchanging the first and third columns on Eq. (|8"Tj) we obtain Eq. ([71)1 . Setting 6 = n 
in this case would lead directly to Eq. (|73|) . The cases with {0n, 0i2, 013} equal to {0,0,0} and {0,0,0} are 
obtained from these through column permutations. 



{011,012,013} -{0,0,0} 





X 




X 


X 


ri = 


X 


, r 2 = 


X 


X 




X X 









(82) 



Performing $1 <H> $2 and exchanging the first and third columns on Eq. ([82)1 we obtain Eq. ([72]) . Setting 6 = n 
in this case would lead directly to Eq. (|75|) . The cases with {0n, 0i2, 013} equal to {0,0,0} and {0,0,0} are 
obtained from these through column permutations. 



{011,012,013} -{0,0,20} 





X X 








r x = 


X X 


, r 2 = 












X 





(83) 



Performing $1 H $2 on Eq. (|83|) we obtain Eq. ([78| . Setting = n in this case would lead directly to the 
special case of = it in Eq. (fTTj) . The cases with {0n, #12, #13} equal to {0,20,0} and {20,0,0} are obtained 
from these through column permutations. 



{0n,0i2,0i3} -{0,0,20} 





X 




X 


r x = 


X 


, r 2 = 


X 




X 




X 



(84) 



Performing $1 «-» $2 and exchanging the first and second columns on Eq. (|84|) we obtain Eq. (|77|) . Setting 
= 7r in this case would lead to Eq. (|73| . after interchanging the second and third columns. The cases with 
{0n,0i2,0i3} equal to {0,20,0}, {0,20,0}, {0,0,20}, {20,0,0}, and {20,0,0} are obtained from these through 
column permutations. 



{011, 012, 013} = {0, 0, 0} 









X 


X 


X 


ri = 


XXX 


, r 2 = 


X 


X 


X 



any 0. 



(85) 



Performing $1 f> $2 on Eq. (f8"5|) we obtain Eq. (fTTj) . Notice that the special case of — ir had already shown 
up in Eq. (|79j) . 



5 Recall that the cases {012,023,031} = {9, — 9,0}) and {012,023,031} = {— 9,0,9} are obtained from this through permutations on 
the rows of the Yukawa matrices. 
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• {011,012,013} = {M, 20} 









X 


X 


r x = 




, r 2 = 


X 


X 




X X 






X 



Performing $ x O < &2 on Eq. ([551) we obtain Eq. ((75)) . Setting = 7r in this case would lead to Eq. (1751) . after 
interchanging the first and third columns. The cases with {0n, 0i2, 013} equal to {0,20,0} and {20,0,0} are 
obtained from these through column permutations. 

For those wishing to check that all possibilities have been considered, we refer to the footnote 6 . 



4- Special cases with 8 — tt 



We continue to explore the cases where each ay can only take the values 0, 0, or —0 (mod 2ir). Certain cases 
are only valid for 8 — tt. For example, consider a\2 — (mod 2tt) and 023 = 031 = (mod 2ir). Taking G [0, 2tt[, 
this can only happen for 9 = tt, due to Eq. (IM1) . This forces us to consider the case {«i2, «23, «3i} = {0,7r, n}. The 
cases {ai2, «23, «3i} = {tt, 0, tt} and {«i2, Q?23, «3i} = {tt, tt, 0} are obtained from this by permuting the rows on the 
respective Yukawa matrices. In this case, 02i = 0n, 03i = 0n + 7T, and 



e 



Cll &12 (713 

011 012 013 

011 + TT 6\2 + 7T 013 + TT 



(87) 



implying that 9\j G {0, 7r}. There are 2 3 such cases, all of which lead to a matrix T where all entries may be non- 
vanishing 7 . We continue to ignore cases which differ only by permutation of the columns. It is easy to see that we 
have already considered all possible structures. Indeed, when all 9\j equal tt, we recover Eq. (|71[) ; when two du equal 
7T, we recover Eq. (|73[) : when only one 9\j equals tt, we recover Eq. (|75l) ; and when no 9\j equals tt, we recover Eq. (|79l) . 



5. Special cases with 8 = 2n/3 

We now turn to the last two cases where each ay can only take the values 0, 0, or —0 (mod 2tt). Due to Eq. (|54l) . we 
can have ai2 = 023 = 031 = ±0 if and only if = 2tt/3. The case ai2 = 023 = a^i = —2tt/3 (or, which is the same, 
47r/3) is obtained by exchanging any two rows of the Yukawa matrices. We choose the case ai2 = 023 = a^i = 2tt/3, 
implying that 02i = 9%i — 2tt/3, 631 = 9%i + 2n/3, and 



e = 



"11 "12 "13 

0H-27T/3 0i2-2tt/3 0i 3 -2tt/3 
0ii+2tt/3 0i2 + 2tt/3 i3 + 2tt/3 



implying that 6\j G {0, 27r/3, — 27r/3}. Recall that 0n = c 12 = C13 
quarks. 

Ignoring cases which differ only by permutation of the columns, we are left with the following structures: 



is excluded because it would lead to massless 



We have also checked that 

• The cases where {0n, $12, $13} equal {0, 0, 0} lead to vanishing quark masses, if 9 7^ tt, and to Eq. H71II . if 9 = tt; 

• The cases where {0 n , 9 12 , 9 13 } equal {0, 29, 29}, {29,0, 29}, and {29,29,0} lead to vanishing quark masses, if 9 ^ tt, and to Eq. J7TJ, 
if 6» = tt; 

• The cases where {0 n , 9 12 , 9 13 } equal {9, 29, 29}, {29,9, 29}, and {29, 29,9} lead to vanishing quark masses, if 9 + tt, and to Eq. ([73}, 
if 6» = tt; 

• The cases where {9\\, 6*12, #13} equal {29, 29, 29} lead to vanishing quark masses, if 9 ^ tt, and to Eq. (1711 . if 9 = 7r. 

Of course, some entry may be zero by accident. The point is that this value is not required by a symmetry of this type and, as such, it 
is not invariant under the renormalization group equations. 



13 



{011,012,0X8} = {0,0, 2tt/3} 





X X 




X 


ri = 


X 












X X 



{011,012,013} = {0,0, -2tt/3} 



{011,012,013} = {0,27r/3,27r/3} 



{011,012,013} = {0,27r/3,-27r/3} 



{0n, 012, 0i 3 } - {0, -27T/3, -2tt/3} 



{0n, 012, 0i 3 } - {2tt/3, 2tt/3, -2tt/3} 



{0n,0i2,0i3} - {2vr/3, -2vr/3, -2tt/3} 



= 2tt/3; 





X X 






ri = 




, r 2 = 


X 




X 




X X 



2tt/3; 





X 




X X 


ri = 


X X 


, r 2 = 










X 



27T/3; 





X 




X 


ri = 


X 


, r 2 = 


X 




X 




X 



= 27t/3; 





X 






ri = 




, r 2 = 


X X 




X X 




X 



2tt/3; 









X X 


ri = 


X X 


, r 2 = 


X 




X 







27T/3; 









X 


ri = 


X 


, r 2 = 


X X 




X X 







= 2tt/3; 



(89) 



(90) 



(91) 



(92) 



(93) 



(94) 



(95) 



Care must be exercised when comparing these matrices with those shown previously. Consider, for example, Eq. (1891) . 
{0n, 012, 0is} = {0, 0, 2tt/3}, with {6»2i, (9 3 i} = {0ii - 2tt/3, 6 n + 2tt/3} = {-2tt/3, 2tt/3}. When might worry about 
Eq. (|T2"j) . where one can also choose {0n, 012, 013} = {0, 0,27r/3}. However, there, {0 2 i,03i} = {0n,0n + 27r/3} = 
{0,2tt/3}. 



D. Yukawa matrices for up- type quarks 



So far, we have only shown the Yukawa matrices for the down-type quarks. We will now show that it is trivial to 
get the Yukawa matrices for the up-type quarks from those for the down-type quarks. Let us start from some specific 
transformation of the left-handed fields, characterized by ai2 and 0:31. From Eqs. (|3"0"j) and (|4"T) we get 02i = 0n — ai2, 
031 = 0ii + c*3i, so that 



e 







fll (712 (713 

711 — a 12 012 — Oi\2 #13 — OL\i 

?11+Q!31 012+O31 013+^31 

011 012 013 

\\ — Ct\2 0_12 — Q!l2 0_13 — a 12 

Ill+Q-ZX 012+^31 013 + a 31 



(96) 
(97) 
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Each entry on the column j of 6 is of the form 9±j + b. We then followed the procedure 

9ij + 6 = (mod 2tt) => entry is in Tx, 

0ij + b = 9 (mod 27r) =► entry is in T 2 . (98) 

Let us call 9\j = 9\j — 8. Then, if 9\j +6 = (9), we find 9\j + b = —9 (0), meaning that this is an entry in A 2 (Ai). 
Thus 

9 X j + b = (fly - 6) + b = -9 (mod 2tt) =4> entry is in A 2 , 

hj + b = (9ij - 9) + b = Q (mod 2tt) =*> entry is in Ai. (99) 

The argument goes both ways, so we can find all cases for the up-type Yukawa matrices A by starting from all cases 
for the down-type Yukawa matrices T and performing the following procedure: 

• 9ij — > 9ij = 9ij — 9; 

• ri — > A 2 ; 

• T 2 — > A x . 

Of course, one can shuffle differently the columns of -{Ti,r 2 } and {A 2 , Ai}, since they live on different right-handed 
spaces. 



E. Counting the number of models 



The only purpose of our parameter counting is to show the enormous amount of cases which have been killed by the 
simple requirements that there be no massless quarks and that the CKM matrix not be block diagonal. As pointed 
out at the end of section HU there are potentially 3 18 = 387.420.489 different models. Notice that this number does 
not include permutations that lead to the same form for the Yukawa matrices. But, it does include permutations 
which, although leading to different forms of the Yukawa matrices, have no impact on the physical observables. These 
same procedure must be followed when we count the number of distinct forms of the Yukawa matrices based on the 
analysis of the previous sections. 

The forms shown in sectionHTfflcorrespond to 6 L x (3+3+3+3+6+6+6+6)„/? x (3+3+3+3+6+6+6+6) pi ? = 7776. 
The sub-indices L, nR, and pR correspond to the permutations of rows, down-type columns, and up-type columns 
(respectively), that lead to the same physics. But, as in the 3 18 possibilities above, the counting has been performed 
so that no two structures look the same. The numbers in (3 + 3 + 3 + 3 + 6 + 6 + 6 + 6) n R correspond to the number 
of possibilities in Eqs. (f5"T]) - (j6"4"]) . respectively. 

To be specific, let us look at Eq. (|57[) . Exchanging the first and second column leaves the form invariant. This 
is counted as one structure. However, exchanging the third and first columns leads to a new structure. So does an 
exchange between the third and second column. There are thus three possibilities. This explains the first "3" in 
(3 + 3 + 3 + 3 + 6 + 6 + 6 + 6) n R- The rest of the counting procedure follows the same lines. 

The forms shown in section UlI C II correspond to l^x (1 + 3 + 3 + l) n R X (1 + 3 + 3 + 1) p r = 64. The forms shown 
in section UTTCl with 9 £ tt c orrespond to 3 L x (1 + 3 + 3 + 6 + 3 + 3) nR x (1 + 3 + 3 + 6 + 3 + 3) pR = 1083. The 
forms shown in section UlI C 21 with 9 = tt correspond to 3l x (3 + 3 + 3)„r x (3 + 3 + 3) p r = 243. The forms shown 
in section ITlI C 31 correspond to 3l x (1 + 3 + 3 + 6 + 3 + 3) nR x (1 + 3 + 3 + 6 + 3 + 3) p r = 1083. Finally, forms shown 
in section UTTCl] correspond to 6 L x (3 + 3 + 3 + 6 + 3 + 3 + 3) nR x (3 + 3 + 3 + 6 + 3 + 3 + 3) pR = 3456. There are 
thus 13.705 distinct surviving possibilities. 

This may seem like a large number, but notice that we have eliminated 387.406.784 a priori conceivable Yukawa 
structures. The simple requirements of quarks with nonzero mass and a CKM matrix which is not block diagonal 
provides a drastic reduction in the number of possibilities. Said otherwise, the huge majority of Yukawa matrices 
consistent with abelian symmetries do not survive simple experimental constraints. We should also point out that any 
two structures which differ only by permutations of the rows (simultaneously in T and A), and/or by permutations 
of the columns of T, and/or by permutations of the columns of A give exactly the same physics. Permutations aside, 
we are left with the 8 + 4 + 9 + 6 + 7 = 34 possibilities for the down- type Yukawa matrices shown in Eqs (|57l) (|64| , 
(|7I ]H [79"] ) . ([5T ]) — ([55 ]) . and dSSJ)-®, with similar structures for the up-type Yukawa matrices. Combining 
appropriately, we get 8x8 + 4x4 + 9x9 + 6x6 + 7x7 = 246 overall models. Those that differ only by $i o $ 2 
will lead to the same physics. Of those, a few can be further excluded because they do not yield any CP violation. 
The possibility of spontaneous CP violation will be addressed in section |V] 
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IV. TWO IMPORTANT RESULTS 



A. Most discrete symmetries have the same impact 

We have considered a symmetry in the scalar sector S — diag{l, e l9 }. Of course, if the lagrangian is invariant under 
S, it is invariant under any power of S. In this way, if 8 = 2tt /n, then the Z n group is generated. If 8 ^ 2-rr/n, then 
one generates a discrete, but infinite, group. For simplicity we will refer to the Z n groups in what follows. 

We now turn to an important result from our previous analysis. We know that choosing 8 — 2ir/3 or 8 — 2-k jh 
leads to the same Higgs potential. Indeed, any 8 ^ 0,7r leads to the same Higgs potential as the continuous U(l) 
Peccei-Quinn symmetry [3[. From this point of view, applying any Z n (n > 3), or even U(l) is the same. With the 
results presented in the previous section, we see that this is no longer the case when the fermions are added. As shown 
here, the symmetry Z3 allows Yukawa structures not allowed for other Z n . Remarkably, all Z n with n > 4 have the 
same impact on the full Lagrangian, even when fermions are introduced. 



B. Most discrete symmetries imply an accidental continuous symmetry 

The notation 8 7^ 2ir/3, tt, Att/3 used in Eqs. (|5T1) - (|64I) means that the form of the matrices shown is the most 
general consistent with values of 8 which differ from 2ir/3, tt, and 4ir/3. But one should notice that the form of the 
matrices shown are left invariant even if 8 = 2tt/3, 7r, 4-7r/3. The point is that, in general, for those special values of 8 
these matrix forms are not the most general consistent with the symmetries. For example, Eq. (j57[) is not the most 
general matrix consistent with 8\\ = 812 = 8, 813 = 28 (mod 2tt) when 8 = tt. That form is shown in Eq. (|75p . 
But one can see that, indeed, Eq. ([57]) is a particular case of Eq. (|75|) . So, Eqs. (|57|) ~ (|64)l are invariant under the 
symmetry for all 0, but they are only the most general forms consistent with the symmetry for those symmetries 
where 8 ^ 27r/3, 7r, 47r/3. The dedicated reader can check this explicitly by comparing these forms with the forms 
presented for the special cases 8 = tt and 8 ~ 2ir/3. 

This has a very important consequence. A matrix form which is invariant under the symmetry for some value of 
8 =/= 2ir/3, 7r, 47r/3 will be invariant under the symmetry for all values of 8, meaning that the Yukawa sector will be 
invariant under t/(l). Since this is also true for the Higgs potential, we conclude that, for the cases in section UlI Bl i) 
Imposing Z2 on the scalars does not imply a larger symmetry, neither in the Higgs sector, nor in the Yukawa sector; 
ii) Imposing Z3 on the scalars implies a continuous symmetry in the Higgs sector, but not in the Yukawa sector; hi) 
Imposing Z n , n > 4 on the scalars implies a continuous symmetry, both in the Higgs sector and in the Yukawa sector. 

The other cases can be analyzed in a similar fashion. For the cases in section MI C II i) Imposing Z2 on the scalars 
implies a continuous symmetry in Yukawa the sector, but not in the Higgs sector; ii) Imposing Z n , n > 3 on the scalars 
implies a continuous symmetry, both in the Higgs sector and in the Yukawa sector. For the cases in sections MI C 21 
and MI C 31 i) Imposing Z2 on the scalars does not imply a larger symmetry, neither in the Higgs sector, nor in the 
Yukawa sector; ii) Imposing Z n , n > 3 on the scalars implies a continuous symmetry, both in the Higgs sector and in 
the Yukawa sector. 



V. SPONTANEOUS CP VIOLATION 



A. Strict two Higgs doublet model 



Let us now look at the possible vacua of a theory with only two Higgs doublets and three fermion generations, and 
their implications for CP violation at the lagrangian level. We are interested in implementations of discrete abelian 
symmetries, like Z n , for which the scalar potential of Eq. (p} can be written as 



V = m^^J^i + m| 2 ^2*2 — mi 2 (<&i®2 + h.c.) 

+ iAxC^i) 2 + iA 2 ($t$ 2 )2 + A3($J*i)(s£* 2 ) + A 4 ($ t 1 $ 2 )(^$i) + iA 5 fe^) 2 + h-c 



(100) 



where all the parameters are real. We have included the soft-breaking parameter mf 2 , taken to be real so that CP 
is not explicitly broken. For a Z2 symmetry - 8 — tt in Eq. (|23[) - the A5 coupling is present in the potential. For 
Z n , n > 3, or indeed any other value for 8 different from or tt, the symmetry sets A5 to zero and the potential is 
indistinguishable from the Peccei-Quinn one [ll| . At the minimum, the scalar fields develop vevs which we take to 
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be given by, without loss of generality 

($ 1 )=vi=u 1 , ($ 2 ) = v 2 = u 2 + iu 3 , (101) 

with all Ui real. A vacuum with 143 7^ may lead to spontaneous CP violation (SCPV) in the scalar sector - however, 
the presence of a phase in the vacuum is no guarantee of SCPV. To verify whether SCPV occurs in the scalar sector, 
we must calculate the basis invariant quantities of ref. [12j, which was done for all possible THDM scalar potentials 
in (l3| . The minimization conditions are given by dV/dui = 0, from which we obtain 

= [m 2 n + \ lU 2 + (A 3 + A 4 )(«2 + uf) + \ 5 (uj - iif)] "i - m 2 12 u 2 (102) 

= [m 22 + X 2 (ul + u 3 ) + (A 3 + A 4 + A 5 )ztf] u 2 — m\ 2 ui (103) 

0= [ml 2 + \ 2 (u 2 2 + ul) + {\ 3 + \ t ~ \ 5 )u 2 ] u 3 . (104) 

From these we see that solutions with u 3 = are always possible. There are several interesting cases: 

• 9 = ir, exact Z 2 symmetry (m 2 2 = 0, A5 7^ 0): from Eqs. (|104[) and f|103[) . any solution with 1*3 7^ automatically 
implies either u\ = or u 2 = 0. Both solutions lead to no SCPV in the scalar sector (see 

• 6 = 7T, softly broken Z 2 symmetry (mf 2) A5 7^ 0): both solutions without SCPV in the scalar sector {u 3 = 0) 
and with SCPV in the scalar sector (u 3 7^ 0) are possible, depending on the values of potential's parameters (l4j . 

• 9 7^ {0,7r}, exact U(l) symmetry (m 2 2 — A 5 = 0): the equations above only determine the sum u 2 + itj, and as 
such the relative phase of the vevs is arbitrary. These vacua lead to no SCPV in the scalar sector [13[ and in 
fact generate an axion. 

• 9 7^ {0,7r}, softly broken U(l) symmetry {m\ 2 7^ 0, A5 = 0): from Eqs. (|104|) and ()103|1 . we see that any solution 
with u 3 7^ leads to u\ = which, considering Eq. (I102[) . also implies u 2 — 0. Thus, no SCPV vacuum can 
occur in this case. Vacua with u 3 = possess no axion. 

The existence of an axion in one of the cases above is easy to understand: as was explained earlier, the imposition 
of a discrete symmetry with 9 7^ {0, 7r} (for instance a Z n symmetry with n > 3) on the scalar potential leads to an 
accidental Peccei-Quinn continuous U(l) symmetry. Any vacuum for which both fields acquire a vev will break that 
symmetry and lead to a zero mass for the pseudoscalar. This corresponds in fact to the appearance of an additional 
Goldstone boson (other than the three usual ones arising from the breaking of the gauge symmetry). Analytically, 
the pseudoscalar mass is given by 

v 2 

m\ = m? 2 + 2A 5 w 2 , (105) 

1*11*2 

with v 2 = u\ + u 2 , for vacua with 1*3 = 8 . From this we see that: the Z 2 potential will never lead to an axion, since 
A5 7^ 0; the exact U(l) symmetry forces this mass to be zero; and the softly broken Z n potential again has no axion, 
as the pseudoscalar mass is directly proportional to the soft breaking parameter. 

The scalar vevs originate the fermion masses, but also have a contribution to CP breaking at the lagrangian level, 
whether they are real or complex. In fact, the Jarlskog invariant, which measures CP violation in the weak interactions, 
is given by [HI 

J = Tr[H u , H d f = &i(m 2 t - m 2 c )(m 2 - m 2 u )(m 2 c - m 2 u )(m 2 - m 2 s )(m 2 - m 2 d ){m 2 s - m 2 )Im (V us V cb V: b V c * s ) , (106) 



where the matrices Hd and H u have been defined in Eqs. (|44[) and (|45j) . In the SM, since no CP breaking can arise 
spontaneously, it is explicitly broken with complex Yukawa couplings. In the THDM we can study models where one 
has demanded that the full lagrangian be CP invariant, such that the matrices Ti and A; will be real, and the only 
possibility of producing a non-zero Jarlskog invariant will be the vevs having a relative phase. Since such a vacuum 
is impossible for the softly broken U(l) scalar potential, we conclude that models with an abelian symmetry (other 
than Z 2 ) and with an explicit CP conservation are ruled out, since for them J would always be zero. Nonetheless, 



In the case of the exact U(l) symmetry an arbitrary phase between the vevs is possible, but it has no effect on the scalar masses 
whatsoever. 
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there is a distinction worth making : the special forms found for the matrices with 9 = 2n/3 (Z3 symmetry), given 
in section 1111 E[ would give a non-zero Jarlskog invariant if a vacuum with a complex phase could be produced; all the 
other Yukawa matrices we have obtained for the cases 9 7^ 7T, 2tt/3 give J — even if a complex vacuum existed. As 
such, the only models allowed are those, like the SM, where CP is explicitly broken by the Yukawa couplings. 

As for the Z2 model, the exact symmetry is also ruled out when CP is explicitly preserved - no phase from the vevs 
can originate J ^ 0, even for the odd case U2 = 0, allowed by Eqs. (j 102[) — (|104|) : in that case there is a phase of ir/2 
in the vaccum, but it has no bearing on J, which gives zero. In the softly broken Z2 model a vacuum with a relative 
phase in the vevs may be obtained and it leads to CP violation, both in the scalar and the Yukawa sectors [lij ■ And as 
before, Z2 models with explicit CP breaking are in principle perfectly viable. We summarize this analysis in Table HI 

TABLE I: Possibilities of CP violation for THDM with abelian symmetries. "Yes" means that the model's parameters can 
generate a non-zero value for the Jarlskog invariant. The "U(l) n models are those for which one has imposed a discrete 
symmetry of the form of Eq. (|23p . with 6 7^ 0, n. 



Model 


Lagrangian with 
explicit CP breaking 


CP-conserving 
Lagrangian 


Exact Z2 


Yes 


No - real vacuum 
or 

vev phase gives J = 


Soft-broken Z2 


Yes 


Yes 


Exact U(l) 


Yes 


No - vacuum 
gives axion 


Soft-broken (7(1) 


Yes 


No - vacuum with phase 
impossible 



A few observations are in order: 

• We have not considered in this analysis the so-called "inert vacua", where either ($1) = or ($ 2 ) = 0, possible 
in the case of exact symmetries (Z2 or U (1)). These give an acceptable J only in the case of explicit CP breaking. 

• The Z3 case is special. Let us again consider the case of explicit CP conservation. Unlike the remaining 
symmetries with 9 7^ ir, 2w/3, a vacuum with complex vevs would give J ^ 0. Such a vacuum is impossible in 
the THDM, but one can conceive (like the authors of @ did) models with two doublets and additional gauge 
singlets, capable of producing the desired form for the vevs jig . 

B. Complex vacua and the Jarlskog invariant 

The vacua of a Z n potential may be easily altered by introducing soft-breaking terms, as discussed in the previous 
section, or by the inclusion of extra singlet scalars. Here we discuss those cases where the introduction of singlet 
scalars implies a relative phase between v% and V2, and we ask whether this provokes the appearance of a phase in 
the CKM matrix when all Yukawa couplings are real 9 . 

To do this we calculated the Jarlskog invariant of Eq. (|106p . assuming a relative phase between v\ and V2 for all 
the 246 models of Yukawa matrices (assumed real) which we have identified. In almost all cases J = 0. The only 
exceptions occur for 9 — tt or 9 — 2ir/3. The results are presented in Table [TTI and Table [Till respectively. These 
tables will be useful for the study of spontaneous CP violation in models with two scalar doublets and various scalar 
singlets, in the presence of abelian symmetries. 



Of course, the inclusion of scalar gauge singlets has no impact on the Yukawa matrices we have found in the previous sections, since 
singlet scalars have no coupling to the fermions. 
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TABLE II: We assume that 9 = n, that all Yukawa entries are real, that the vevs have a relative complex phase, and we 
calculate J. The down-type Yukawas were chosen according to the equations along the first line, and the up-type Yukawas 
were chosen according to the equations along the first column. We denote the entries where J — 0, all others allow for J / 0, 
depending on the values of the parameters. 



Equations 










for Yukawa 




(|73|) 


(|75|) 


d79|) 


matrices 










CD 











m 










m 










(79) 












TABLE III: We assume that 9 — 2tt/3, that all Yukawa entries are real, that the vevs have a relative complex phase, and we 
calculate J. The down-type Yukawas were chosen according to the equations along the first line, and the up-type Yukawas 
were chosen according to the equations along the first column. We denote the entries where J = 0, all others allow for J / 0, 
depending on the values of the parameters. 



Equations 
for Yukawa 
matrices 


(|89|) 


(|90|) 


m 


(|92| 


(193} 


(|94| 


(|95|) 


(ESI 


















in 


















(pTJ) 


















(|92]) 
















((93]) 


















ami 


















(E0> 



















VI. NATURAL SUPPRESSION OF FLAVOUR CHANGING NEUTRAL SCALAR INTERACTION 

Measurements in the mixing of neutral mesons (such as K — K, Bd — Bd, etc.) lead to tight constraints on 
flavour changing neutral scalar interactions (FCNSI). The discrete symmetry Z 2 was introduced in the scalar sector 
by Glashow and Weinberg @ and, independently, by Paschos Q, precisely to preclude such FCNSI. But there are 
several other options to curtail FCNSI. For example, one may invoke large scalar masses, or introduce approximate 
flavour symmetries [8j. Perhaps more interestingly, one may relate the FCNSI with the CKM matrix. In a very nice 
article, Branco, Grimus, and Lavoura (BGL) used discrete abclian symmetries in order to construct one such THDM 
@, following earlier work by Lavoura [||. The BGL model corresponds to the use of our Eq. (|78[) for the up-type 
Yukawa matrices and of our Eq. (|71[) for the down-type Yukawa matrices. 

One may now ask the question: is there any other implementation of abelian symmetries which leads to a relation 
between FCNSI and the CKM matrix? Although we have all possible implementations of abelian symmetries, the 
question is difficult to answer analytically because it involves diagonalizing the mass matrices. Indeed, the quark mass 
basis is obtained with the basis transformation 

d-L = V d L n L , 
dn, = VdR riR, 

UL = V U L PL, 

UR = V U R PR, (107) 

where we have used = (tt,l,pl) t . The unitary matrices VdL, VdR, V u l, and V u r are chosen such that 

d\d.g{m d ,m s ,m b } = D d = V d L [v\Ti + v 2 T 2 ] Vj R , 

diag{TO„, m c ,m t } = D u = V uL [v{ Ai + v* 2 A 2 ] V^ R . (108) 
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The CKM matrix is V = V uL Vj L . The matrices controlling the FCNSI are 

N d = VdL h*r! - <r 2 ] vJ R , 

N u = V uL [v 2 A, - Ul A a ] VI R . (109) 

Botella, Branco, and Rebelo [13] have proposed a method to identify BGL-type implementations while sidestepping 
the diagonalization procedure. They start from the relation [B[ 

N d = ^-D d - —v dL r 2 vJ R , (110) 

obtained by combining Eqs. (| 108[) and (|109[) . and using v 2 = \vi\ 2 + \v 2 \ 2 . Based on this they propose the following 
sufficient conditions for BGL implementation: i) Ai + v 2 A 2 is block diagonal; and ii) there exists a matrix P such 
that ha) PT 2 — kT 2 (for some number k), and iib) PTi = 0. As they stress, the condition can be applied with an 
up- type/down- type quark interchange. 

We start by noticing that Eqs. (|T08| and (|T09l) can also be combined into 

„2 



N d = — —D d + —VdL Ti Vj R , (111) 
v 2 v 2 



implying that an equally good sufficient conditions for BGL implementation is: i) u*Ai + v 2 A 2 is block diagonal; 
and ii) there exists a matrix P such that ha) PI\ = ATi (for some k), and iib) Pr 2 = 0. Again, the condition can be 
applied with an up- type/down- type quark interchange. The new condition is just a 3>i -f-> $ 2 transformation of the 
previous condition, useful to us when looking for all possible BGL-type implementations. 

Since we have tabled all possible matrices, we are able to see that only Eq. ([78)) can lead to a block diagonal 
v{A\ + v 2 A 2 for the up-type quarks. We must now check all compatible down-type Yukawa matrices, namely, 
Eqs. dnl), (|72]l. (|74|) . (|76l) . dZZ]), and see whether they satisfy condition ii) 10 . We have checked that only for Eq. ([7T|) 
can one find a matrix P consistent with the constraints ii). 

This gives a unique character to the work of Branco, Grimus, and Lavoura |9j. They have developed the only 
possible implementation of a relation between FCNSI and the CKM matrix which uses abelian symmetries and is 
consistent with the sufficient conditions above. There are only two caveats. First, we have only checked the sufficient 
conditions developed by Ref. [Kj and extended here. A priori, one can entertain the possible existence of cases which 
do not satisfy the sufficient conditions presented, but where the FCNSI are indeed related to the CKM matrix. In the 
cases where we could perform the analysis analytically, we have found no such case. Second, in some cases condition 
ii) is violated because it leads to constraints on the non-zero matrix elements of the Yukawa matrices. It could be that 
some non-abelian group might lead to further zeros on the Yukawa matrices, thus evading the problem. Although 
possible, such a case would be difficult to construct because more zeros in the Yukawa matrices will, more often than 
not, lead to massless quarks or to a block-diagonal CKM matrix. 

In light of our analysis, that a BGL Q case was found by inspection in the THDM is truly remarkable. 



VII. CONCLUSIONS 



We have studied the restrictions on the Yukawa matrices imposed by discrete abelian symmetries acting on the 
scalar and fermion sectors of the THDM. Using known experimental constraints, we have reduced the number of 
possible cases from 3 18 to 246. Ignoring row and column permutations, we are left with 34 types of down-type 
Yukawa matrices (and the same for up- type quarks), which we table explicitly. 

We have found that imposing a symmetry Z n (n > 4) on the scalars always leads to an accidental U(l) symmetry; 
that applying a Z3 symmetry on the scalars leads to an accidental U(l) symmetry in the scalar sector but not 
necessarily in the fermion sector; and that applying a Z 2 symmetry on the scalars does not lead to an accidental U(l) 
symmetry in either sector. 

We show that only Z 2 with soft breaking in the scalar sector enables spontaneous CP violation. We also show that 
the proposal of Branco, Grimus and Lavoura [9( is unique, in our context, and conjecture that this uniqueness might 
hold even when non-abelian symmetries are considered in the THDM. 



The possibility that both the up-type and down-type Yukawa matrices are given by Eq. (1 78 1) is excluded, since it would lead to a 
block-diagonal CKM matrix. 
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Finally, we stress that our results have a very wide applicability in model building because all discrete non-abelian 
groups have a Z n subgroup, for some value of n. For a given non-abelian group, pick one of its Z n subgroups and 
diagonalize its generator. Applying that generator as a symmetry of the lagrangian, one falls into one of the 34 Yukawa 
matrices we have shown explicitly. The action of further generators (which, of course, need not be diagonalizable in 
the same basis) will, in general, lead to further constraints on the Yukawa matrices. Given the low number of entries 
in many of our Yukawa matrices, and the likelyhood of further constraints setting them to zero, the action of further 
generators will often lead to matrices inconsistent with experimental constraints. 
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